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By means of new general variational method we report a direct solution for the quintic self-focusing 
nonlinearity and cubic-quintic ID Gross Pitaeskii equation (GPE) in a harmonic confined potential. 
We explore the influence of the 3D transversal motion generating a quintic nonlinear term on the 
ideal ID pure cigar-like shape model for the attractive and repulsive atom-atom interaction in Bose 
Einstein condensates (BEC). Also, we offer a closed analytical expression for the evaluation of the 
error produced when solely the cubic nonlinear GPE is considered for the description of ID BEC. 
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I. INTRODUCTION 

Nowadays one and quasi-one dimensional Bose- 
Einstein condensates (BEC) are common experimen- 
tal procedures^. The transition from 3D to ID sys- 
tem was invoked long time ago^. In general, the 3D 
Gross-Pitaesvkii equation (GPE) cannot be factorized 
into transverse and longitudinal motions, nevertheless, 
under certain parameter regions we can assert that the 
BEC follows a ID behavior (for a detailed discussion see 
Ref. 0). In the case of the harmonic trapping potential 
and considering that the atoms are tightly confined in 
two transverse directions, a transition to the quasi-lD de- 
scription is possible. Starting with the standard 3D GPE, 
employing the adiabatic approximation and using the 
anzat wavefunction &(x, r;t) = exp(— inot/h)$(x)x(r',t), 
we can derive an effective ID GPE, which describes the 
physical characteristics of the cigar-like shape conden- 
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where € R is the chemical potential, m is the atomic 
mass, u) is the longitudinal harmonic oscillator frequen- 
cies, and gm € R, g € R are the effective ID nonlin- 
ear self-interaction coefficients. These two coefficients 
depend on the total number N of particles in the con- 
densate, the transverse harmonic oscillator frequency uj r 
and the scattering length a s (a s > or a s < for attrac- 
tive or repulsive interatomic interaction, respectively) by 
the relations g 1D = 2a s Nhuj T and g = 6\n(4/3)g^ D /huj^., 
where we have chosen for the stationary state 4> the nor- 
malization condition J dx |$| 2 = 1. 

H 

Equation P is a cubic-quintic nonlinear Schodinger 
equation (NLSE) with real coefficients. The presence 
of the — g\<&\ term in ([T]) is due to the deviation from 
one dimension on the longitudinal condensate dynamics, 
i.e. a residual three dimensionality on an effective one- 
dimensional GPE. In the case of a homogeneous medium, 
i.e. assuming that uj is zero, the cubic-quintic NLSE was 



widely used to describe the physical process of an optical 
medium with a nonlinear polarization including suscep- 
tibilities up to fifth order—. Also, the soliton solutions 
have been extensively studies in Refs. [5081. 

As we have mentioned above, Eq. (JlJ is the cigar-like 
shape approach from the 3D GPE. An important issue is 
the range of validity of (JXJ) , which is directly linked to the 
existence and stability of set of ground states solutions of 
the 3D NLSE. It is well known that for any value of a s > 
the 3D GPE does not collapse^. However, for attractive 
interatomic interactions, the solution is dynamically sta- 
ble if and only if a s is within the ranged 
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< 0.627, 



(2) 



with a±_ = yJh/mui T . Hence, the validity of the cigar-like 
shape approach represented by Eq. (JlJ is also restricted 
to the constrain 

Rescaling to dimensionless variables l = y/h/mu, 
£ = x/l , A = 2g 1D /(l hu), fi = 2/i /(Mi V>(f) = 
®(&o)/Vh, e = 31n(4/3)cj/w r , Eq. §Q can be cast as 



M 4 V = /^- (3) 



The main task of this paper is the implementation of a 
more general variational mathematical approach to solve 
Eq. ([3|). Based on this result, we provide approximate 
solutions for the order parameter, the chemical poten- 
tial and minimal energy for the quintic and cubic-quintic 
GPEs. The paper is organized as follows. First, in Sec. 

II we present the bases of our formalism, i.e., we present 
exact formulae for the energy and the chemical potential 
as functions of relevant parameters of (j3)). By consider- 
ing a trial function for the ground state, we derive in Sec. 

III a representation for the energy, £ , app (A), and chemical 
potential, /i a pp(A). Section IV is devoted to the applica- 
tion of our results to get explicit approximate solutions 
for the quintic and also cubic-quintic NLSEs. An estima- 
tion of the error due to the influence of the interaction 
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between the axial and radial degrees of freedom on the 
ID cigar-shape model is presented both graphically and 
analytically as function of the self-interaction parameter 
A and the coefficient e leading the quintic nonlinear term. 



II. CUBIC-QUINTIC NONLINEAR 
GROSS-PITAEVSKII EQUATION 

In the following we will consider a more general non- 
linear Gross-Pitaevskii equation 



-ft+e^ + aX^l^-bX'^ri;^^, ^€1, (4) 

where a > and b are real constants. 

Let V= UjeH\R)\ /£ 2 |V(£)| 2 ^<+°°} be the 
Hilbert space endowed^ with the norm 
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and the corresponding inner product 



W#)v = 



where 
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\mfdti and \m)\\l 



(8) 

are the usual norms of the standard Banach spaces L 4 (]R) 
and L 6 (R) . Taking into account the relation between the 
energy and the chemical potential, i.e. SE\[tfj\)/8ip\ — 
fi^x, it follows that 



Stpx ' d\ 
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dX 
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Hence, from Eqs. ([7]) and (|9]) we obtain the useful formula 



dExbh 
dX 
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Thus, the minimum energy is given by 



A A 
£min(A) = l + | f \\ip s \\ 4 4 ds-j [ S \\lP s \f 6 ds. (10) 











Moreover, from Eq. ^ is straightforward that the chem- 
ical potential can be written as 



In V we define the energy functional 
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We denote by Gx the set of ground states of Eq. ((H), 
i.e., the set of functions of V that minimize the energy 
functional Ex[ifi] under the condition 



(6) 



Notice that in the case of attractive interaction where 
b > 0, it is possible to show (by applying the Gagliardo- 
Nirenberg inequalities) that the set of ground states Gx 
is nonempty if the condition |A| < 2/y^jbf is satisfied. 
Hence, for V>a € Gx we obtain 



dExbP. 



SEx&x] d^x\^a 4 26A 6 
dX \^ x ->^x/ + 2^- — ^> 

(7) 



Mmin(A) = E min {X) + — Wi'xWi g- WxWq ■ 



or equivalcntly 

Atmin(A) = 1 + 7 

_2b 
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(11) 



It is important to remark that Eqs. (fTUj) and (jlip arc 
exact under the condition of knowing the ground state 
0a G Gx and therefore, independent of the method or 
approach we employe to get the solution of the Eq. <j3j> - 



III. APPROXIMATE FORMULAE 



It is possible to showiS that any solution of Eq. dU 
belonging to the Hilbert space V has the asymptotic be- 
havior exp(— t£ 2 ), with r > 0, as |£| — > oo. So, to evalu- 
ate the minimal energy -E m i n (A), we can consider for the 
ground state the trial function 
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(- 
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exp(-r£ 2 ). 



(12) 



Using the function (fl"2")l and evaluating the energy func- 
tional ([5]), we obtain the algebraic expression 



E x [il> T ] = (1 - eA> 2 + ^2 



aA 



(13) 



with e = 2&/(3v3V) and cr = ^/r. In the case e < the 
Eq. (1131) presents a global minimum in (0, +oo) for any a 
and A € M., but for e > a global minimum is guaranteed 
if |A| < = v/W3/V2&. 

Let cr(A) = y r(A) be the minimizer of Eq. (fT3")) . in 
this way the ground state solution in C/^ can be searched 
considering the function <p\{£,) := VV(A)(£)- By using 
Eqs. (fTUj) and (fTTj) we obtain the approximate energy, 
E app {\), and chemical potential, jU a pp(A), namely 



E app {X) — 1 + 



20? 



cr(s)ds-2e / scr 2 (s)ds (14) 



and 



Mapp(A) = 1 + 



2y^ 



Acr(A) + / cr(s)ds 
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AV 2 (A) + / sa z (s)ds 
o 



(15) 



It becomes clear that |A| < 2/y/b < Ag = y / 37rv / 3/ '\/26, 
ensuring that the set C/a is nonempty and that the func- 
tions E app {\) and /^ a pp(A) can be good approximations 
of £min(A) and /u m in(A), respectively. 



IV. APPLICATIONS AND DISCUSSION OF 
THE RESULTS 

Explicit formulae (fTU)) , (ITT1) and the approximate re- 
lationships (IT41 . (fT5")) are among the main results of our 
work. Nevertheless, more explicit expressions for the or- 
der parameter <px , the energy and the chemical potential 
as a function of the atom-atom interaction term would 
be desirable. To do so, let <7 m i n = a min (a,e;X) be the 
global minimizer of Eq. flT3j). Thus, Eqs. ([14]) and (fl"5|) 
with a = (Jmin allow to obtain the dependence E app 
and fj, a pp on the relevant physical parameters a, e and 
A. Depending on the values of A and the sign of the pa- 
rameter e, we can distinguish several phases linked to 
both type of interaction strengths, i.e. i) pure attractive 
(A < 0, e > 0), ii) pure repulsive (A > 0, e < 0), and hi) a 
mixture (A < 0, e < or A > 0, e < 0). 



More precisely, the value of T min (a min (A) := ^r min ) 
that minimizes the function r i — > E\[ip T ] fulfil the equa- 
tion 



(1 - eA > 4 + 
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Notice that, for e < or e > under the condition |A| 
< 1/ y/e, the function (fl~3|) is strictly convex and coercive 
on the interval (0, +oo), and in consequence Eq. (ITB1) has 
a unique solution, while for |A| > 1/y/e the function (fT3|) 
is not bounded from below. 



A. Quintic NLSE 

Firstly and for sake of comparison, we will consider 
the quintic NLSE in a harmonic potential^. Choosing 
a = 0, from (fll))) we have 



2 



1 



(17) 



2(l-eA 2 ) 1/2 ' 
Inserting $F7§ into Eqs. CH]) and flU) with a = 0, we 
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FIG. 1: (Color online) /i^s — A map for the parameter values 
e — ±0.06 and ±0.012. According to the signs of A and e 
the character of the fj. v 5 , as approximated solution of Eq. Q 
with a = 0, can be mapped into fourth zones: A > 0, e < 
-pure repulsive interaction; A < 0, e > -pure attractive 
interaction; A > 0, e > and A < 0, e < -mixture region. 

obtain for the approximate energy E v n and chemical po- 
tential /j, v 5 the expressions 



E v , = (1 - eA 



2\l/2 



and 



AV = - e ^ 2 ) 



2\V2 



eA 2 



(l-eA 2 ) 1/2 ' 



(18) 



(19) 
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Figure Q] displays the fi^s — A map diagram for several 
values of the parameter e. It can be seen that the re- 
duced chemical potential ^i ¥ ,5=2/Lto/(Sw) shows a strong 
dispersion as a function of A, moreover and following the 
symmetry properties of Eq. (|19l) , we observe that in- 
creases (decreases) for the pure repulsive phase, A > 0, 
e < (pure attractive phase, A < 0, e > 0), while in 
the mixture region an opposite behavior is reached with 
respect to the pair of values (A, e). 

Following Eqs. (fT2"j) and (fTTj) we obtain for the wave- 
function <Ptp5 \ the expression 



<p<p b ,x(£) 



i 



?r(l - eA 2 ) 



1/2 



1/4 



exp 



2(l-eA 1 



,1/2 



(20) 

valid for 1 > A 2 e. The above obtained wave function 
exhibits different behavior depending on the sign of e 
and independent of the type interaction (attractive with 
A < or repulsive for A > 0). The function becomes 
effectively less confined for e < 0, i.e. <p v s t \ is delocalized 
and its maximum decreases, while for e > the function 
if »s a(0 g e ts more localized and the maximum increases 
as the nonlinear potential A 2 e increases. 



B. Cubic-quintic NLSE 
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FIG. 2: (Color online) Reduced chemical potential /i v 3_ v 5 as 
a function of the dimensionless self-interaction parameter A, 
calculated from Eqs. (I22[l - (|27[ l and several values of coefficient 
e. Inset: A - £ map showing the validity ranges of Eqs. (|21|l- 
(|27|) and the stability region of 3D GPE given by — A < As = 
2.857/(e) 1/2 (orange) and -A < X 3D = 2.33/ (e) 1/2 (gray 
lines), respectively. 



Although one can solve Eq. (TTol) numerically and to 
obtain from Eqs. (fT4|) and (fl5)) the energy E (p 3_ (p 5(X,b) 
and the chemical potential jtx ¥ ,2_ ip 5(A, 6), it will be very 
useful to report explicit compact approximate solution of 
the cubic-quintic nonlinear ID GPE (QJ. Searching the 
solution of Eq. (fl6|) with a = 1 and b = e > as a Taylor 



series on A, cr 



we get 
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Under the condition |e| << 1 and by substituting 
Eq. (j2"Tj) in (fT4"|) and (TT51) we have, for A small enough 



E v s (A) 



\/2A 1 1 

Ay^TT IT 64 
1 1 
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V2 



1277^256 



A 



-A 



V2 



32tt 2 256 2 13 tt 2 0F2 6 



A 



(24) 
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The terms AE ip a_ ip 5 (A, e) and A^j, ip a_ ip & (A, e) tell us how 
much error is made by neglecting the potential — eA 2 
in the ID cigar-like shape model. These quantities are 
expressed by 



E 3 



AV 3 



AE, 



(A,e), 



^a(A) + Afj, (p a_ v> s(X,e) 



(22) 



(23) 



where E lp 3(X) and /^(A) correspond to the energy and 
the chemical potential, respectively, for the cubic NLSE 
and are given, up to the 5th order, byiS 



AE lf 3_^(X,e) 



V3 
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-A + 



1 
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eA, (26) 
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Figure is devoted to the calculated chemical potential 
(j, lp 3_ lp s using Eqs. (|23j) . (|25|) , and ( |2T|) as a function of A 
for e = 0, 0.05, and 0.1. First, by comparing Figs. (JTJ and 
© we see the strong qualitative difference between the 
two types equations here considered, quintic and cubic- 
quintic NLSEs. Notice that this difference is remarkable 
even at A « 0. From the Fig. © we can assert that 
the cigar-like shape approximation retaining term up to 
1 0| 2 is a good approach solely for the repulsive case. The 
chemical potential (and also, the energy) for A > and 
lj << uj r is almost independent of e = 3 ln(4/3)o;/w r , 
while a not negligible contribution is reached to /i ip 3_ ip 5 
if the atom-atom interaction is attractive, even for very 
small value of e. In the former case the term — eA 2 \(j>\ 4 is 
responsible for the strong obtained dispersion compared 
to the e = limit. In order to understand from the phys- 
ical point of view the behavior of the chemical potential 
on A, we define the effective potential 




FIG. 3: (Color online) The potential U e ff given by Eq. J28]) 
for A = 5 (solid lines) and A = —5 (dot-solid lines). In the cal- 
culation we employed the values of e = (green), 0.05(bluc) 
and 0.1 (black). The chemical potential /J, lfi 3_ lp 5 for each e is 
indicated by flat lines. 



and repulsive interatomic (A = 5) interactions. Also, in 
the figure is represented the values of /j lp 3_ (p 5 for e = 0, 
0.05, and 0.1. It becomes clear that for a given A > 0, 
the variation of the function U e ff(£; A, e) with respect to 
the parameter e is negligible. Hence, the corresponding 
chemical potentials /x„3_„5 (A, e) ~ ^3_^b(A,0). Thus, 
the nonlinear potential A \ ip\ 2 for A > becomes a very 
good description to study the physical properties of a ID 
cigar-shape BEC under the condition that the transverse 
harmonic oscillator frequency u r is much larger than the 
perpendicular frequency trap u. This result is in agree- 
ment with the conclusions of Ref. [l2|. As one can see 
from Fig. ([2]) the error Afj, lp 3_ lp 5 (A > 0, e) ranges between 
2-7 % with respect to the value of ^,3 (A). Now, if we 
consider the opposite case, i.e. an attractive interaction, 
the behavior of [7 e //(£;A, e) presents a strong localized 
potential and, in correspondence, the chemical potential 
will change drastically as the parameter s increases. This 
characteristic of the effective potential (|2"51) for A < de- 
termines clearly that the term A \^\ is not enough for a 
correct description of ID cigar-like shape BECs. In this 
case, the residual three dimensionality term — eA 2 \ip\ , 
for the effective ID GPE under a harmonic trap, strongly 
modify the corresponding chemical potential. Small vari- 
ation of the strength eA 2 leads to strong change of fi v 3 



value and hence the ground state energy £L 



as well. 



This peculiar behavior is related to the orbital stability 
of the nonlinear Schrodinger equation (j4|) . As mentioned 
above, the 3D GPE presents a set of ground states which 
is orbitally stable for any value of the self-repulsive inter- 
action, while for the attractive interparticle interaction 
regime (a s < 0) the solution does not collapse if and only 
if the condition Na s > — aj_0.627 is fulfilled. The for- 
mer criterion can be rewritten as —A < X^d = 2.33/v / £ 
and it should be compared with condition of existence 
and validity, for —A < As = 2.857/y^, of the obtained 
solutions Eqs. p2]) -(p7 ]l . In the inset Fig. (0 for the at- 
tractive interaction, the intersection of these two sets, 
A3u(e) n As(e), is represented by a shaded region in the 
A — e diagram. 

In conclusion, a new variational approach is presented, 
which allows to construct for the cubic-quintic GPE 
closed analytical expressions for the order parameter, the 
minimal energy, and the corresponding chemical. Focus- 
ing on the compact analytical expressions, we report the 
contribution of the quintic term and the systematic er- 
ror of the residual 3D contribution to the ID cigar-shape 
model for both attractive and repulsive interaction. By 
the calculations the obtained chemical potential solution 
highlight the strong dependence on the sign of strength 
A and the values of quintic self-interaction parameter s. 



^//=£ 2 + A|^a| 2 -£A 2 |^| 4 , (28) 

where the order parameter ip\(£) := VV(A)(£) has 
been substituted by the trial function (fT2"|) with 
a (p 3_ {p 5 (A; e) := V /T ip 3_ ip s . Figure © shows the poten- 
tial U e ff for both considered cases, attractive (A = —5) 
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